Electric Control of the Exciton Fine Structure in Non-Parabolic Quantum Dots 
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We show that the non-parabolic confinement potential is responsible for the non-monotonic be- 
havior and sign change of the exciton fine-structure splitting (FSS) in optically active self-assembled 
quantum dots. This insight is important for the theoretical understanding and practical control by 
electric fields of the quantum state of the emitted light from a biexciton cascade recombination 
process. We find that a hard-wall (box) confinement potential leads to a FSS that is in better 
agreement with experimentally measured FSS than a harmonic potential. We then show that a 
finite applied electric field can be used to remove the FSS entirely, thus allowing for the creation 
of maximally entangled photons, being vital to the growing field of quantum communication and 
quantum key distribution. 
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Introduction. Entangled photons, a non-classical state 
of light, are an indispendible part of proposed and im- 
plemented protocols for optical quantum communication 
and quantum key distribution pQ. Parametric down- 
conversion (PDC) is a well-studied and established way 
of creating entangled photons but suffers from two major 
drawbacks: since on the order of 10 10 pump photons are 
required per created photon pair, the process is rather in- 
efficient. Secondly, the creation time is highly stochastic. 
Most protocols, however, require an on-demand source of 
fixed photon number, leading to the search for alterna- 
tives to PDC. One of the most promising candidates is 
the biexciton cascade recombination from a single semi- 
conductor quantum dot (QD) 2J considered in this paper. 

The biexciton is the bound state of two electrons and 
two holes in a semiconductor. Upon radiative recombi- 
nation, one electron and one hole are annihilated by the 
emission of a single photon, leaving a single exciton. As 
a next step in the cascade, the exciton recombines and 
emits a second photon. Because of the Pauli exclusion 
principle the ground state of the biexciton is a singlet, 
whereas there are two possbile exciton states, character- 
ized by the underlying p x or p y orbitals of the hole. Thus 
the emitted light can have either horizontal, \H), or ver- 
tical, \V), polarization when emitted from the biexciton 
and from the intermediate exciton. If the two exciton 
states are energetically indistinguishable, the two recom- 
bination paths are equivivalent with respect to the fre- 
quency of the emitted light, and the two photons are 
entangled in polarization, being described by the state 
(\HH) + \VV))/V2. The process is schematically shown 
in Fig. [lj 

For real QDs, however, the intermediate exciton states 
are typically not degenerate, but split by an energy S 
known as the exciton fine-structure splitting (FSS). Since 
the two recombination paths now have energetically dif- 
ferent intermediate levels, a simple frequency measure- 
ment could reveal the "which-way" information, and thus 
the entanglement is, at least partially, lost [3]. Under- 
standing of the FSS is essential for any model consider- 
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Figure 1. The bixeciton recombination cascade. For degen- 
erate intermediate exicton levels (S <C F, where T is the 
linewidth of the emitted light and 5 the FSS) the two decay 
paths are indistinguishable and the recombination can be de- 
scribed by | XX) |0) -> (|X) |H) + |Y) \V))/s/2 -> |0) (|HH) + 
\W))/y/2 where we use the notation | electronic) photonic) 
and \XX) is the biexciton state, \X) (\Y)) denotes the ex- 
citon state with horizontal (vertical) polarization, \H) (\V)) 
the photonic state with horizontal (vertical) polarization and 
|0) is the electronic ground state or photon vacuum. 



ing the entanglement of the light produced by the biex- 
citon cascade recombination and there are several differ- 
ent sources contributing to the FSS, which have been the 
topic of a number of experimental and theoretical stud- 
ies in recent years. These include intrinsic effects caused 
by the underlying crystal orientation [IHZ], piezoelectric 
effects [8 and the QD geometry 9 11 . There is a funda- 
mental physical difference between the intrinsic splitting 
and geometry dependent splitting. As pointed out by 
Singh and Bester [7], for the intrinsic case, the FSS de- 
creases for smaller dots, as opposed to the geometry de- 
pendent effect, where a smaller dot gives rise to a larger 
FSS. We claim that these two statements do not stand in 
conflict with each other, since they describe two different 
mechanisms. The intrinsic splitting is caused by the ac- 
tual shape of the underlying atomic orbitals whereas the 
geometry dependent one is mediated by the electron-hole 
overlap. This means that the intrinsic FSS depends on 
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the number of atomic orbitals, which increases for larger 
dots, in contrast to the geometry dependent FSS which 
scales with the electron-hole overlap that decreases with 
increasing dot size. In this paper, we presume that the 
crystal orientation has already been chosen to eliminate 
the intrinsic FSS and focus on the geometry dependent 
splitting caused by the long-range electron-hole exchange 
interaction and determined by the in-plane shape of the 
QD. We especially investigate the possibility of restoring 
degeneracy by virtue of applying an in-plane electric field. 
The dependence of the FSS on an in-plane electric field 
has been experimentally studied and it was confirmed 
that the FSS can be removed P21 OS] but a complete 
theoretical understanding is still missing. Earlier work 
using a harmonic dot confinement potential find qualita- 
tively different results [TU]. Here we use a more realistic 
confinement potential and compare the obtained FSS to 
the one from harmonic confinement and to experimental 
results. Our theory explains the observed FSS in InGaAs 
QDs including the sign change of the FSS when a lateral 
electric field was applied [13] ■ In addition to the appli- 
cation of lateral electric fields, there are other methods 
of reducing the FSS by strain [5], vertical electric fields 
[T7HT9] or a magnetic field [20], as well as proposals to 
recover the entanglement in the presence of a finite FSS, 
such as spectral filtering [TU [JS] , time reordering [T5] or 
embedding the QD in an optical cavity [2"T] . 

Theoretical model. We consider a quantum dot com- 
posed of one cubic semiconductor surrounded by a ma- 
terial composed of another cubic semiconductor. The 
electronic structure is characterized by bands of which we 
consider the conduction band states, labelled m s = ±1/2, 
and the valence band states, labelled rrij = ±3/2 for 
heavy holes and rrij = ±1/2 for light holes. Because of 
the large distance in energy to these valence band states 
we do not consider the split-off band. The surrounding 
material is required to have a conduction(valence) band 
of higher(lower) energy at the T-point than the dot ma- 
terial. We start from a general Hamiltonian, describing 
two particles in a semiconductor, 

H Q = H 1 {v 1 ) + H 1 {v 2 ) + V{r 1 -r 2 ), (1) 

where H\ (r) is the single particle Hamiltonian, 



of two-particle systems an anti-symmetric wave function 
is needed and thus we form the two-particle basis 
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Inserting this into Eq. ([I]) , we find a k dependent Hamil- 
tonian containing an interband k • p term. Treated as a 
perturbation up to second order, these terms alone give 
rise to the effective mass approximation. As shown by 
Pikus and Bir [22 , including the Coulomb exchange be- 
tween particles leads to more correction terms such as 
the band diagonal first order Hartree correction, 

^kii'k^skjtks = ^ao'^tt' V r k' 1 _k 1 5ki+k a ,k' 1 +ki> ( 5 ) 

with Vk the Fourier transform of the Coulomb potential 
as well as the third order exchange term, H a , with ele- 
ments 
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where E g is the band gap energy, a and (3 run 
over the spatial coordinates, and (ri, r 2 |ik 2 ski) = 
( / 3 s k i(r2)</5tk 2 (ri)- With the exchange interaction known 
for basis vectors, the two-electron wave functions can be 
expressed as 

* st (ri,r 2 ) = 22 c skl ,tk 2 ^k 1 tk 2 (ri,r2) 
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We choose 



Hl(r) = j~ + lattice (r) 
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with the underlying periodic lattic potential, lattice (r), 
and V(ri — r 2 ) = e 2 / Aneoe r \ri — r 2 |, the Coulomb poten- 
tial with the relative dielectric constant e r . In the k • p 
method a basis for the electron wave functions is formed 
from the Bloch waves as 



¥> s k(r) = e lkr u s (r), 



(3) 



where s is a band index and u s (r) = u St k=o(r) is the 
Bloch function at the T-point of band s. For the study 



where 



^t(ri,r 2 ) = ip 3 (ri)ipt(r 2 ), (8) 



V> s (r) = 5> sk e k ' r Ws (r) = F s (r)u s (r), (9) 



and F s (r) is the envelope function which is the solution 
without corrections ^ and ^ to the equation 
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where m* is the effective mass of band s and V s (r) is a 
band dependent potential varying on a mesoscopic scale 
describing the confining structure such as quantum dot. 
The new correction terms, ([5| and Q can now be ex- 
pressed in terms of envelope functions and are equal to 
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A similar expression was presented by Kadantsev and 
Hawrylak [10] as well as by Tong and Wu [11]. If we 
let t' ,t represent valence band electrons we need to re- 
verse the order and apply the time-reversal operator 
when going over from the electron-electron picture to the 
electron-hole picture. This amounts to 
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where is the time-reversal operator [22] . 

A suitable potential, V s (r), has to be choosen to prop- 
erly describe the system under consideration, in our case 
the quantum dot. We discuss different choices in the next 



section. For a given V s (r), we numerically solve Eq. (10) 
to find the envelope functions. This provides us with 
the single particle states \s — m s ) and \t = rrij) for elec- 
trons and holes from which exciton product states are 
formed as \m s ,mj) = \m s ) ® |m ? -). Using these envelope 
functions, the integrals Eqs. (fTTJ and (12) are calculated 



numerically and the following exciton eigenvalue problem 
is formulated in the basis of \m s , m 9 ) as 



(H Q + H C + H a ) \X) =E X \X) 



(14) 



with H° being the excitation and confinement energies 
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in which t s , s£ TUj} are the energies from Eq. ( 10 1 



This eigenvalue problem is solved numerically and two of 
the eigenvectors, \X X ) and \X y ) arc identified as the ones 
having maximal projections on |cr+) + |<t_) and |<r + ) — 
|(7_) respectively, where \a±) = |=Fl/2, ±3/2) for heavy 
excitons and |er±) = | ±1/2, ±1/2) for light excitons. The 
fine structure splitting is now calculated as 



E x — E„ 



(16) 



with exciton energies E x ,E y taken from Eq. (14). 

Results. We now discuss two choices of confinement 
potentials, V s (r). The objective is to describe the con- 
finement of electrons and holes to a given nanostructure. 
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Figure 2. The FSS 5 calculated for an harmonically confined 
InGaAs dot with charecteristic lengths 30 x 20 x 7 nm 3 and 
various directions of the electric field. Regardless of the di- 
rection of the field, the FSS is monotonically decreasing and 
vanishes only asymptotically. 



In our case we are interested in a quantum dot with di- 
mensions l x xl y xl z . Typically we let l x ^ l y which causes 
a FSS. Physically the confinement comes from a change 
in the underlying composition such as going from InGaAs 
to GaAs. It is necessary to take a band-dependent po- 
tential into account since electrons and holes are usually 
subject to different band offsets. Further, we also in- 
clude an electric field which is represented by — q s ~E ■ r 
in the Hamiltonian, where q s is the charge of a particle 
in band s, i.e. — e(+e) for conduction (valence) bands. 
Here, we restrict ourselves to heavy, bright excitons, i.e. 
|Tl/2,±3/2). 

Harmonic potential. A simple model for the QD is the 
harmonic confining potential 



^ r 2 



,E-r. 



(17) 



The solutions are harmonic oscillator wave functions 
with characteristic lengths l sa — yjh/m*ui sa defining the 
spread of the wave function and lf a — q s E a /m* s u> 2 a , the 
electric displacement. The FSS in this case as a func- 
tion of the electric field is plotted in Fig. [2] We observe 
that the FSS is decreasing with increasing field, but does 
not change sign. Experiments, however, show another 
picture including a non-monotonic behaviour as well as 
a change of sign [13]. This is because in a hamonically 
confined dot, the electron and hole can be arbirarily sep- 
arated only due to the electric field. This neglects actual 
"hard wall" confinement of the dot, which is not affected 
by the electric field, and this suggests another model is 
necessary to properly understand the experimentally ob- 
served FSS. 

Hard-wall confinement. To incorporate the effects of 
a physical confinement we a consider rectangular box of 
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Figure 3. The FSS 8 calculated for Ir10.2Gao.sAs/GaAs dot of 
dimensions 55x50x7 nm 3 with a hard- wall confinement in the 
presence of an electric field applied in various directions. The 
FSS is sensitive to the direction in which the field is applied: 
In the ir-direction the FSS decreases and changes sign at a 
critical field strength where 5 = (here, at E x ~ 35kV/m). 
When the field is applied in the y-direction, the FSS increases 
instead. 



dimensions l x x l y x l z which has a potential step at the 
boundary, i.e. 
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where AE S is the band offset and xa(t) is the charac- 
teristic function of the set A. As can be seen in Fig. [3] 
the FSS now exhibits more structure and depends more 
drastically on the direction of the applied field. A com- 
parison between the envelope functions found for the 
two different cases reveals that the harmonically confined 
particles are displaced but not deformed by the electric 
field, whereas for the case of a hard-wall confinement, the 
shape of the envelope function is modified as well. Fig. [4] 
shows the effect for the hole envelope function. When an 
electric field is applied, the harmonically confined parti- 
cles move outside physical boundaries of the dot, whereas 
the particles confined by a step potential remain inside 
the box. By inspecting Eq. (12 1 we also note that ex- 



change integral depends not only on the electron-hole 
overlap but also on the curvature of the envelope wave 
function. The curvature, in turn, depends on E in the 
case of the hard- wall confinement, but not for the har- 
monic confinement. 

Conclusions. We theoretically investigated the exciton 
FSS for quantum dots of cubic semiconductor materials 
and its dependence on dot geometry and applied electric 
field. We found that the choice of confinement poten- 
tial in the model is of great importance and we noted 
that different confinement potentials can lead to quali- 
tatively different results. A model with a harmonic con- 
finement potential cannot capture the experimentally ob- 
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Figure 4. Hole envelope functions with, F^.(r), and with- 
out, F mj (r), applied electric field E x in x-direction, for (a) 
harmonic and (b) hard-wall confinement. For the case of har- 
monic confinement, the wave functions are translated by Z^. 
and are no longer inside the physical region of the dot when 
a field E x is applied but retain their shape. For the hard- 
wall potential, the main effect of the electric field is to deform 
the wave functions; this is accompanied by a relatively small 
shift that leaves the particles inside the dot. The deformation 
affects the second derivative of (r) which determines the 



FSS, see Eq. (12 I 



served features including a sign change of the FSS under 
the application of an electric field. Using a more realistic 
potential step (hard wall) barrier, we find a more complex 
relation between FSS and field and predict the possibility 
of a complete suppression of the FSS, as observed exper- 
iments [13] . We trace the additional FSS structure back 
to the fact that electron and hole wave functions are not 
only displaced but also deformed by the hard-wall poten- 
tial. This deformation influences the FSS via the second 
derivative of the envelope function. The suppression of 
the FSS by means of an electric field allows the creation 
of entangled photons without additional post processing 
which is of interest to the field of quantum communica- 
tion. Open questions include the effects of an applied 
vertical field as well as the influence of light holes. 
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